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The atom optics of Bose-Einstein condensates containing a vortex of circulation one is discussed. 
We first analyze in detail the reflection of such a condensate falling on an atomic mirror. In a 
second part, we consider a rotating condensate in the case of attractive interactions. We show that 
for sufficiently large nonlinearity the rotational symmetry of the rotating condensate is broken. 
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I. INTRODUCTION 

During the last years, the experimental achievement 
of Bose-Einstein condensation (BEC) ||] has raised a 
large interest, and numerous theoretical and experimen- 
tal efforts have been devoted to the analysis of its intrigu- 
ing properties. This interest is partially motivated by the 
fact that the BEC is a macroscopic coherent matter wave, 
with evident applications in the already well-developed 
field of atom optics j3|. In this sense, the dynamics of 
BECs interacting with atom optical elements, such as 
e^g. optical or magnetic mirrors or wave guiding 

H , has been recently investigated. 

Additionally, the BEC presents remarkable properties 
due to its superfluidity , in particular the quantization 
of vortex circulation In this sense, studies of vortices 
in trapped condensates have attracted a growing atten- 
tion H |h], A first aim of this paper is the analysis of 
the atom-optics of a BEC that contains a vortex. In par- 
ticular wc study a BEC with a single vortex of circulation 
one when falling to and being reflected from an atomic 
mirror. We show that the coherent dynamics exhibits in 
such a situation self-interference effects, restoration of a 
vortex of opposite circulation after reflection, splitting of 
the BEC into two parts at the top of the bounce, and for- 
mation of additional vortices in low density regions. This 
analysis may serve as a way to investigate the properties 
of rotating BECs. 

During the recent years, the development of the 
Feshbach-resonance technique has allowed it to change 
the strength and even the sign of the interparticle in- 
teractions |l2|, This technique consists of employ- 
ing magnetic fields to excite resonances between atomic 
and molecular states, and in this way strongly influence 
the value of the s-wave scattering length a s [Q. By 
using this novel technique, particularly remarkable ex- 
periments have been performed, including the Bose-nova 
experiments at JIL A [[l3| , or the recent creation of bright 
solitons (lj, |l6| . The same technique can be employed to 
analyze the stability and eventual collapse of a BEC con- 
taining one or more vortices. This constitutes the second 
part of this paper. We show that a larger stability, sim- 
ilar to that discussed for the case of surface modes [[l7| , 
can be expected in the case in which the rotational sym- 
metry is kept, due to the centrifugal barrier imposed by 
the vortex. However, we show that the rotational sym- 



metry is actually broken when the attractive nonlinearity 
is adiabatically increased to a sufficiently large value. For 
that case, we show that eventually a piecewise collapse is 
produced for nonlinearities comparable with the critical 
ones for a non rotating condensate. 

The paper is organized as follows. In section |l| we 
discuss the reflection of a rotating BEC from a mirror. 
In section III we study the collapse of a rotating BEC, 
discussing the stability and the possibility of breaking 
of the rotational symmetry. Finally, in section IV we 
present some conclusions. 



II. REFLECTION OF A ROTATING BEC 

In the following, we consider a gas of N bosons placed 
in a harmonic trap with pancake symmetry of frequen- 
cies tu±_ -C u> z . For sufficiently strong axial confine- 
ment, such that n -C hu> z , with p, the chemical po- 
tential, the axial dynamics can be considered as effec- 
tively frozen. Hence, the wave-function can be written as 
ip3T>(%i V, z ) = i>(%i y)'0o( z ): where the transversal profile 
ijjQ (z) is that of the ground-state of the axial harmonic 
oscillator. For sufficiently low temperatures, the BEC 
dynamics is then well described by the two-dimensional 
Gross-Pitaevskii equation (GPE): 
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with g 2D = VirH^D, where g 3D = ^-Na s , and m the 
atomic mass. Demanding by constraint, that the non- 
interacting BEC has a vortex of charge n in its center, 
the lowest-energy state in an axially symmetric trap is 
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The reflection of a non-rotating BEC has been exper- 
imentally accomplished both from optical Q and mag- 
netic j^] mirrors. In the following we analyze such a re- 
flection for the case of a rotating condensate. 

We consider in the following a hard mirror (at y = 0), 
corresponding to a sufficiently large and abrupt potential 
barrier. In general the mirror has a finite razor, which 
by means of holographic techniques in the case of optical 
mirrors can be limited to few laser wavelengths, typically 
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of the order of few microns. This length scale is smaller 
than those involved in our case (see Figs. |l| and ||), and 
therefore the softness of the mirror is not expected to in- 
troduce any significant modification in the discussed ef- 
fects. In this case, the mirror can be easily simulated by 
forcing the wave function to be antisymmetric on the axis 
perpendicular to the mirror: ip(y) = —ip(—y). Instead of 
one BEC, which is reflected by a potential, we may then 
consider two condensates: the "real" BEC, and its anti- 
symmetric counterpart on the other side of the mirror, 
the "ghost" BEC. During the reflection both condensates 
just pass through each other, the ghost becoming the ob- 
served reflected BEC. 

In our numerical simulations we have set the trap fre- 
quencies to uj± = 2ir x 20 Hz and u> z = 2tt x 800 Hz, which 
guarantees the two-dimensional character of the dynam- 
ics for the cases considered. By means of imaginary-time 
evolution |D| we create a BEC with a vortex in its cen- 
ter, which is initially located at a height y = h = 78.5 
/im over the mirror. At t = 0, the trap in the xy plane is 
switched-off (ui± = 0). Notice that the axial trapping po- 
tential remains switched-on, and therefore the system is 
maintained two-dimensional. The condensate falls down 
in the gravity field towards a hard mirror at y = 0, where 
it is bounced upwards again. We discuss the results for 
two cases: i) a very low number of atoms (gzo = 0), 
and ii) N ~ 18600 23 Na atoms, which corresponds to a 

coupling parameter g 2 o = 345.6^-. In Figs. |l| and || we 
depict the evolution of the falling BEC at different times 
for the i) and ii) case respectively. The condensate is 
dropped at t = ms, and the reflection occurs at t ~ 4 
ms (the mirror is along the bottom border of the boxes), 
whereas the upper turning point is reached at t ~ 8 ms. 

In the pictures taken at t = 4 ms an interference pat- 
tern can be observed of the falling parts of the condensate 
and the already reflected parts. It is not affected by the 
nonlinearity and shows up both in linear and nonlinear 
simulations. The interference lines are very dense and 
therefore the pattern is not very clearly visible in the fig- 
ures. Shortly after t = 6 ms the BEC is focused to a 
small strip by the gravitational cavity. At the top of the 
bounce the condensates splits into two parts (see expla- 
nation below) with a region of zero density separating 
them. In the interacting case the core size is smaller, 
and consequently the low-density region which separates 
both parts is narrower, as observed in Fig. ^|. In addition, 
as commented below, in the low-density line we have ob- 
served in the interacting case the creation of additional 
vortices p(J. 

The focusing of the condensate at one point in its 
bouncing has classical origins. At the reflection grav- 
ity changes from an accelerating force to an decelerat- 
ing force. Due to the fact that the lower parts of the 
condensate are reflected earlier than the upper parts, all 
atoms gain a small velocity boost towards the center of 
the BEC. Indeed, if the velocity spreading Av of the con- 
densate satisfies Av <C \/2gh, then one can consider the 
condensate particles initially with v = 0. In that case if 



the initial condensate spatial width Ay <C h, one obtains 
that after the first bounce the particles are focused at 
y ~ 3/i/4 at a time t = 3^h/2g + 0((Ay/h) 2 ). Similar 
arguments can be applied for larger velocity spreadings. 
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FIG. 1: Reflection of a rotating BEC for g2o = 0. The gas is 
released at t = (top-left figure). The mirror is placed at the 
bottom line of the figures. 



From the density profiles it is difficult to draw con- 
clusions whether a vortex remains in the BEC after the 
reflection. Much more information can be obtained b 
looking at the phase of the wave function. In Fig. 
a plane wave is added to the wave function, ip{f) — » 

ip{r) + \^{r)\e ik - p . In this way the phase becomes vis- 
ible and a vortex can be detected as a fork in the image 
|]l9| . In order to be able to observe the phase in the di- 
lute outer regions, a constant value ( 4j of the maximum 
density) is added to the density. 

In the first picture of Fig. || the initially created (right- 
spinning) vortex shows up as a fork which is opened to 
the bottom. Shortly after the reflection (t = 5 ms) the 
fork is opened to the top, which means that the vortex 
has changed its direction of spin during the reflection 
and now rotates to the left. Surprisingly, three vortices 
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t = 7 ms 



t = 8 ms 



FIG. 3: A closer look at the phase of the BEC in Fig. |. 
Vortices are visible as fork patterns: Right-spinning vortices 
show up as forks opened to the bottom, while left-spinning 
vortices can be observed as forks opened to the top. 
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FIG. 2: Same as Fig. § but with g 2D = 345.6^ 



are visible in the last picture (t = 8 ms). Two of them 
must have been created during the squeezed state be- 
tween t — 6 ms and t = 7 ms. Even more interest- 
ing is the fact that the total vorticity is not conserved 
during the squeezing: While at t — 5 ms the BEC con- 
tains one left-spinning vortex, at t = 8 ms this vortex is 
joined by two new vortices which are both right-spinning. 
This vortex instability does not appear in simulations, 
in which the nonlinear interactions have been switched 
off, and therefore requires nonlinear interactions between 
the atoms. Similar instabilities of vortices have been ob- 



served by Garcia- Ripoll et al. 1 20 1 for vortices in trapped 
condensates. 

A rather simple model can be employed to account 
for the basic features of the condensate reflection in the 
absence of interactions. In that case, the dynamics is 
provided by a Schrodinger equation with a potential 
V{y) = mg\y\ (employing the ghost-BEC picture dis- 
cussed above), where g is the gravitational acceleration. 
It is possible to obtain a semi classical approximation by 
splitting the process into three steps: i) expansion while 
falling towards the mirror, ii) reflection, iii) propagation 



off the mirror. 

Ignoring the gravitational effects, which will be in- 
cluded later, and evolving from an original wave-function 
ipo(S), the evolution reduces to a simple expansion in 



FT" 



exp(-?ifc 2 t/2m)FT[V'o] 



free space: ip\{x ,t) 

where FT and FT^ 1 denote the direct and inverse Fourier 
transform, respectively. The main effect of the reflec- 
tion step ii) is a velocity boost depending on the po- 
sition in the condensate. The lower parts of the con- 
densate are reflected earlier than the upper parts, and 
this time difference is essential because at the reflection 
gravitation changes from an accelerating force to a de- 
celerating force. We can classically calculate the result- 
ing velocity boost, Au(y) = — 2y/ti, where t\ = y / 2h/g 
is the reflection time of the BEC center of mass, and 
y is the relative position in the frame of the conden- 
sate. After applying this velocity boost to the con- 
densate, ip2{x) = ipi(x , t\) exp(— imy 2 /Mi), we evolve 
in the stage iii) again in free space until a final time 
ti after the reflection of the center of mass: ipf(x) = 

FT -1 (exp(-ihk 2 t 2 /2m)FT(ilj2)y Fig. | shows the evo- 
lution of the BEC wave-function obtained by means of 
the previously discussed approach, which agrees very well 
with the numerical results shown in Fig. fy. We expect 
that this should be the case as long as the velocity spread- 
ing of the condensate Av <C \[2gh. 

The previous simplified picture allows for an easy un- 
derstanding of the breaking up of the condensate into two 
pieces observed in our numerical simulations. In its way 
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up after the reflection the condensate possesses a velocity 
field which results from the addition of the vortex veloc- 
ity field and the linear velocity boost resulting from the 
gravitational field, v = (u>y/r 2 ,—uix/r' 2 — 2y/tl), where 
u> is the angular velocity. As a consequence the compres- 
sion is not symmetrical around the vertical coordinate of 
the center of mass, y — 0, but on the contrary the left 
part of the cloud concentrates in the upper half-plane, 
whereas the right part does it in the lower one. As a 
result the condensate splits into two pieces. 

t = ms t = 4 ms 



o 



66 I 

-12.5 







t = 6 ms 




t = 8 ms 




46.4 i I 66 I i 

-12.5 12.5 -12.5 . . 

FIG. 4: FT results for the same case as in Fig. ED. 



III. IMPLOSIONS OF ROTATING BECS 

While in the previous sections the interparticle inter- 
actions were assumed to be repulsive, we now consider 
clouds of atoms with attractive interactions, i.e. a s < 0. 
For small enough trapped condensates, the dispersion in- 
duced by the zero-point oscillation of the trap can pre- 
vent the nonlinear focusing provided by the attractive 
mean-field, and in this way a metastable BEC can be 
created However, for sufficiently large condensates 
in two- and three-dimensional trapping geometries, the 
gas is unstable against collapses. Actually, in physical 
situations, the formation of a singularity is avoided due 
to the appearance of two- and three-body losses at large 
densities |2l], [2^, ^3) . It is the aim of this section to dis- 
cuss the physics of rotating attractive BECs, including 
stability and the possibility of symmetry breaking. 

As a first approach, we consider a vortex at the center 
of a two-dimensional rotationally-symmetric BEC, and 
assume that the rotational symmetry tjj{f , t) — x(r, i)e Ic ^ 
is kept for any value of the nonlinearity. As we show at 



the end of this section this assumption actually breaks 
down when the value of |<?2d| (<?2D = — |<?2d|) is adiabat- 
ically increased to a sufficiently large value. 

Our analysis has been performed by means of numer- 
ical simulations of the GPE ([!]). We start with a rotat- 
ing condensate in the absence of interactions, (72D = 0, 
which is provided by (Q) with n — 1, in the presence of 
a central vortex. Then, the coupling parameter |<?2d| is 
adiabatically increased in order to establish the criterion 
for condensate stability [glj . Physically the value of <?2D 
can be modified cither by increasing the number of atoms 
for a fixed negative value of the scattering length, or by 
reducing a2D for a fixed number of atoms by means of 
Feshbach resonances. Note that if, on the contrary, 52D is 
changed rapidly, additional structures as those discussed 
in Ref. |2^] could eventually appear. 

We would like to note at this point, that as long as 
the BEC density is not very large, the GPE (|l|) should 
describe very well the condensate dynamics. However, as 
discussed above, at large densities two- and three-body 
losses do play an important role in preventing the for- 
mation of a singularity. This could be described by in- 
cluding in the GPE the corresponding cubic (two-body) 
and quintic (three-body) damping terms as discussed in 
Rcfs.pl], g3|. It is not the purpose of this paper to an- 
alyze the physics after the collapse occurs, and therefore 
we constrain ourselves to the use of the GPE (0). 

In the later stages of the collapse the trapping potential 
becomes negligible compared to the mean-field energy 
and to the kinetic energy. Therefore (even for initially 
non-symmetric confinements) the later stages of the col- 
lapse are characterized by a single length scale l(t) which 
has a universal scaling law l(t) oc y/t~—to, where to is the 
time at which the BEC collapses into a singular point 
p4j . By using this self-similar scaling on the numerical 
grid it is possible to perform the simulation up to the very 
fi nal st ages of the collapse. Fig. || shows the behavior of 
yj (r 2 ) as a function of the coupling parameter (72D for the 
cases with and without a vortex at the trap center. This 
condensate width presents the expected scaling at the fi- 
nal stages of the collapse. As observed in Fig. |[ if the 
rotational symmetry is preserved for any value of |<?2d| 
the critical values of the coupling constant, <? c , for which 
the collapse occurs, significantly differ between the case 
without and with a vortex. In particular, the centrifugal 
force due to the vortex would amount for a larger stabil- 
ity of the condensate, and consequently for a larger value 
of |<7 C |. We have numerically obtained that in the pres- 
ence of vortex the value of g c is roughly four times larger 
than that in the absence of it: g c (n = 0) = —5.85— and 

q c (n = 1) = —24.15—. A similar stabilization mecha- 
nism was discussed in Ref. |17| for the case of attractive 
BEC in the presence of surface modes. 

An estimation of the critical coupling can be obtained 
by considering a Gaussian ansatz 
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FIG. 5: Mean square radius V< r 2 > for a rotating (solid) 
and non rotating (dotted) two-dimensional BEC if the rota- 
tional symmetry is preserved, as a function of p2D, which is 
adiabatically decreased as linear function of time. 



and minimizing the corresponding energy functional pf 
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In this way we obtain q c (n = 0) = — 27r— Rj —6.28—, 

and g c (n = 1) = -32e~ 2 7r|^ » -13.6^. The Gaussian 
ansatz provides a good agreement with the numerical re- 
sults for the case without vortex, but underestimates the 
stability for the case of n = 1. This difference may be 
explained by the departure of the wave- function from the 
Gaussian ansatz, as we have observed in our simulations 
for the case of coupling constants close to g c . 

In our previous numerical simulations, as well as in 
the analytical estimations, we have assumed that the ro- 
tational symmetry is preserved when the value of g 2 D is 
adiabatically changed. In the final part of this section we 
shall show that fully two-dimensional simulations clearly 
indicate that such assumption must be revised, since the 
rotational symmetry is actually broken in the presence 
of any slight disturbance, when 52D is adiabatically in- 
creased. 

In order to analyze the effects of small perturbations of 
the rotating BEC, we consider a simple model, in which a 
Gaussian-like condensate with a centered vortex acquires 
a slight asymmetry in the density distribution around the 
line x(= r cos (f) = 



ip(r, 4>) 
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where a controls the width of the wave function, and d 
is a small parameter that controls the asymmetry of the 
condensate. Due to the fact that the sign of d can be 
chosen arbitrarily, dE/dd = at d — 0, where E is de- 
fined in Eq. ffl) . The second derivative at d = provides 



information about the stability of the condensate: 
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The value for a, which minimizes the energy at d = is 
given by 



dE 



d=0 



1 87T(7 4 (jj 2 — 8tt — g 2 n 
4 er 3 7r 



0. 



Substituting this expression into Eq. ^ we obtain 

1 6g 2 D 
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d=0 



For > (72D > — |tt, this second derivative is positive and 
the condensate is expected to be stable. For 52D < — §tt 

however, ^r|d=o is negative. Therefore, it becomes clear 
even from this very simple model, that beyond a given 
value of |<72d| the rotational symmetry breaks down. 
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FIG. 6: Two-dimensional simulation of a BEC in a vortex 
state when the value of <?2D is adiabatically decreased. The 
lengths are in units of ^Jh/mu, and Q2u in units of 



Figure |sj shows snapshots taken from a two- 
dimensional simulation. A small asymmetry ( corre- 
sponding to d = 0.01 in Eq. (|) ) is applied to the 
initial state. Starting at <?2D = 0, the coupling parameter 
is then adiabatically decreased until the condensate col- 
lapses. Once a critical value of the coupling parameter is 
exceeded (|<72d| — 7.8), the asymmetry of the condensate 
grows and the rotational symmetry breaks. In particu- 
lar, the density concentrates into two oppositely-located 
peaks. If |<72d| is further increased, the width of the peaks 
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FIG. 7: The distributions y/< (r • ei) 2 > (solid line) and 

a/< (r • e 2 ) 2 > (dotted line) of a BEC in a vortex state, 
where e*i and &2 are the unit vectors pointing into the di- 
rection parallel or perpendicular to the symmetry breaking, 
respectively, as indicated in Fig. |^. The vertical dotted line 
indicates the onset of instability. 



decreases and eventually the peaks collapse (\g c \ ~ 9.7). 
Note that the value of g c is by a factor of about 1.7 
larger than that expected for a non-rotating BEC, but 
since the BEC actually splits into two parts, the critical 
nonlinearity in each peak is approximately the same as 
for the BEC without vortex. 

The emergence of the symmetry breaking can be ob- 
served in more detail in Fig. which shows the value 
of y/ < (f ■ ei) 2 > as a function of g2B, where e% (ea) 
denotes the direction parallel (perpendicular) to the axis 

T 2 

of the symmetry breaking. For <?2D > — 7-8^- the ro- 
tational symmetry holds and the condensate is perfectly 
described by those calculations which employ such a sym- 

h 2 fi 2 

metry. For —7.8— > </2D > — 9.7— the rotational sym- 
metry is broken, and the condensate splits into two parts 
orbiting around the trap center. At this stage the larger 

is |<?2d| the narrower are the peaks of the distribution. 

?i 2 

Finally, for 52D < — 9.7— the condensate is unstable and 
collapses into two points. 



We would like to stress that a very similar picture 
has been observed due to numerical inaccuracy even in 
the absence of any initially imposed perturbation. This 
strongly suggests that any sort of physical noise, e.g. 
thermal one, will lead to a breaking of the rotational 
symmetry when the collapse is adiabatically approached. 



IV. CONCLUSIONS 

In this paper we have analyzed some relevant phenom- 
ena occurring in the atom-optical physics of rotating con- 
densates. In a first part we have analyzed the reflection 
of a rotating BEC from an atomic mirror. We have ob- 
served that the vortex is preserved after the reflection, 
both in the presence and in the absence of interactions, 
and that its rotation is inverted. The combination of 
gravitational effects and the mirror produce the squeez- 
ing of the BEC cloud after its bouncing, and leads even- 
tually to the breaking of the vortex and the formation 
of a notch in the density. In the interacting case, we 
have observed the formation of additional vortices in the 
region of the notch. 

In the second part of this paper, we have considered 
the case of a rotating BEC in the presence of attrac- 
tive interactions. We have shown that if the rotational 
symmetry is conserved, the centrifugal barrier induced 
by the vortex amounts for a larger stability of the sys- 
tem. However, more careful calculations show that, if the 
value of the nonlinearity is adiabatically reduced begin- 
ning by the noninteracting case (either by increasing the 
number of condensed atoms or by modifying the value of 
the scattering length via Feshbach resonances), the ro- 
tational symmetry is actually broken, and eventually a 
piecewise collapse is produced. Such a symmetry break- 
ing should occur for any initial slight noise, i.e. even in 
the absence of any imposed initial external perturbation. 
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